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Nucleus-nucleus collisions : what we have learned from the heavy ion 
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In the first part, I give a brief description of the quark-gluon plasma search at CERN and of some experimental 
results. In the second part, I review a dynamical model of nucleus-nucleus interactions and propose a physical 
interpretation of those results. 
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1. QUARK-GLUON PLASMA SEARCH 

Heavy ion collisions is a very active field in par- 
ticle physics. At present the maximal energy is 
reached at CERN where lead ions are accelera- 
ted to an energy of 160 GeV/c per nucleon. The 
next machine will be the RHIC collider, expected 
to operate in 1999, in which gold or lead ions 
will reach a c. of m. energy of t/s — 200 GeV. 
One decade later an energy of ^/s = 5.5 TeV 
will be reached at CERN-LHC. The aim of these 
programs it to produce and detect a new state 
of matter called quark matter (QM) or quark- 
gluon plasma (QGP). Since ordinary matter is 
composed of quark and gluons it is natural to 
think that by heating and/or compressing it, a 
new phase of deconfined quarks and gluons can 
be reached. This idea is quite old and was first 
explored in the framework of the bag model. Ho- 
wever, it only gained full credibility when lattice 
calculations showed that such a phase transition 
does occur in statistical QCD. 



old. Schematically, the situation is as follows. 
For T < T c one heats the system and the tem- 
perature increases whereas the ratio e/T 4 , where 
e is the energy density, remains constant. At 
T = T c one can heat the system but the tem- 
perature does not increase (only the latent heat 
does). Finally, for T > T c the temperature in- 
creases again and e/T 4 remains constant. This 
new constant is much larger than the correspon- 
ding one for T < T c . This is due to the larger 
number of degrees of freedom of a QGP as com- 
pared to a pion gas. Using the value of the new 
constant and the value of T c one finds that the 
phase transition takes place when the density of 
the system is an order of magnitude larger than 
that of ordinary nuclear matter. It is also found 
that there exists a restoration of chiral symmetry 
which takes place at the same temperature T c . 

When quarks are introduced, the calculations 
are more complicated and the order of the phase 
transition depends on the number of light quark 
flavors. 



1.1. Lattice QCD results 

In the case of a pure Yang-Mills theory (i.e. 
QCD without quarks) the phase transition is of 
first order and takes place at a critical tempera- 
ture T c ~ 150 MeV. In the Big Bang theory the 
same transition occurred, in the opposite direc- 
tion (i.e. from QGP to ordinary nuclear matter), 
when the universe was only a few microseconds 
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1.2. QGP production 

Although QGP could be produced in the case 
of neutron stars we are interested here in its pro- 
duction in heavy ion collisions. The first remark 
is that, in this case, we are a priori very far 
from the conditions of thermodynamical equili- 
brium assumed in the lattice calculations in sta- 
tistical QCD. However, as far as the energy den- 
sity is concerned, it is plausible that values of the 
order required by theory (i.e. 1 GeV/fm 3 ) can 



be reached in central Pb-Pb collisions at present 
CERN energies. In the Bjorken scenario, the e- 
nergy density is given by 

_ [dN/dy} y ^ <E T > m 

7T Rf, b T 

where dN/dy is the number of secondaries per 
unit rapidity, < Et > is the average transverse 
energy per particle, and To is the plasma forma- 
tion time. This time is very poorly known and 
there is no clear consensus among experts in the 
field. Using To = 1 fm together with the experi- 
mental values of the other quantities in ([!]), one 
gets e~ 2 -=-3 GeV/fm 3 . 



slightly above the critical temperature and, se- 
cond, c and c being very rear will combine with 
light quarks to produce open charmed mesons 
and therefore, once melted, the J/ip will not be 
formed again. 

1.4. Experimental situation 

The two phenomena described above have been 
observed experimentally in a very clear way. In- 
deed it has been found that the values of the ra- 
tios 

K A A 3 



1.3. QGP signals 

Even if QGP is produced, the phenomenon is 
a transient one, and one has to find experimental 
signals of its formation. Although some signals 
regarding the bulk properties of the production 
have been proposed in the past (i.e. behavior of 
< pt > versus dn/dy), it seems at present that 
the best signals correspond to rear events or to 
electromagnetic processes. I will only discuss here 
two of these signals : strangeness enhancement, 
and J/ip suppression. 

Due to chiral symmetry restoration it is ex- 
pected that the ratio of strange over non-strange 
particles will increase due to a reduction of the 
production threshold in dense hadronic matter. 

The suppression of J/ip production, proposed 
by Matsui and Satz, is a consequence of color 
screening in a plasma (Debye screening) - similar 
to the screening of electric charge in an ordinary 
plasma. This screening produces a modification 
of the potential of the type 

V{r) = V {r)exp(-r/r D (T)) (2) 

where rjj is the Debye radius - which decreases 
when T increases. When r£>{T) is smaller than 
the binding radius, the corresponding bound state 
cannot exist. Of course, at sufficiently high T, 
there will be a transient melting of all hadrons. 
The interest of the J/ip (as well as other char- 
monium and bottonium states) is twofold. First, 
lattice QCD calculations show that the melting 
of J/ip occurs at temperatures T ~ 1.2 Tc only 



increase between pp or pA and central AB colli- 
sions. The NA35 collaboration claims that there 
is no increase between pp and pA, but errors are 
still large. 

J/ip suppression has been observed by mea- 
suring the ratio J/ip over Drell-Yan (DY). It is 
found that this ratio decreases continuously from 
pp to pA to AB collisions. Furthermore, for a 
given AB system, the ratio decreases with in- 
creasing centrality. This effect is very strong in 
Pb-Pb collisions (see Section 2.6 and Fig. 8). 

1.5. Physical interpretation 

There are claims in the literature that 
strangeness enhancement and particularly the en- 
hancement of multi-strange antibaryons indicates 
the formation of QGP. Nevertheless, it has also 
been possible to explain these data in the frame- 
work of string models - at the cost of introducing 
some new ingredients such as string fusion and 
final state rescattering of secondaries (see section 

2.5) . Likewise, alternative explanations of the ob- 
served J/ip suppression have been proposed (see 

2.6) . 

In the following I will briefly review an indepen- 
dent string model : the dual parton model (DPM) 
and will discuss to what extent the above signals 
can be understood in this framework, what are 
the required modifications of the model and what 
these heavy ion experiments teach us regarding 
the mechanism of multiparticle production. 



2. STRING MODELS AND QGP SI- 
GNALS 



The dual parton model (DPM) is a dynami- 
cal model aimed at describing the mechanism of 
multiparticle production in hadronic and nuclear 
collisions. A very similar model, the quark gluon 
string model (QGSM) has been introduced by 
the ITEP group. These models are based on the 
large- N expansion of non-perturbative QCD. A 
comprehensive review can be found in ref. [Q. 
A short description of the physical basis of the 
model and some of its applications to hadronic 
collisions can be found in [^j. Here I will concen- 
trate on the generalization of the model to nu- 
clear collisions, and, in particular, on the eventual 
modifications of the model required to explain the 
CERN heavy ion data. 

2.1. The model 

In pp collisions the main contribution consists 
in a color exchange producing two strings of type 
diquark quark (see Fig. 1). The hadronic spec- 
tra of each string is obtained by a convolution 
of the momentum distribution and fragmentation 
functions of the constituents at the string ends. 
Due to energy-momentum conservation, the lon- 
gitudinal momentum fractions of the constituents 
at the string ends in each hemisphere add up 
to unity. There are also more complicated con- 
figurations, corresponding to higher order terms 
in the large- N expansion, involving 4, 6, . . . etc 
strings, with sea quarks and antiquarks at their 
ends. These configurations correspond to multi- 
ple inelastic scattering in an S-matrix approach. 
In pp collisions, these higher order configurations 
give small contributions at moderate energies, but 
their contribution becomes increasingly impor- 
tant with increasing energy. In proton-nucleus 
and nucleus- nucleus collisions, these higher-order 
contributions are enhanced by trivial combina- 
torial factors (proportional to A 1 / 3 in pA colli- 
sions) . More precisely, in pA collisions the domi- 
nant graph, involving only one struck nucleon of 
A and A — 1 spectator nucleons, is the same two 
string diagram of Fig. 1. A double inelastic colli- 
sion, involving two struck nucleons of A and A — 2 
spectators, corresponds to a four string graph as 



Figure 1. Dominant two string component in 
nucleon-nucleon scattering. 



shown in Fig. 2, and so on. The weights of these 
contributions can not be computed in QCD and 
are taken from the Glauber-Gribov approach. 

One assumption of DPM, which is not moti- 
vated by the large-A^ expansion, is that strings 
are independent and particles emitted in different 
strings are uncorrelated. This assumption works 
well in pp, pA and peripheral AB collisions. Ho- 
wever, we shall see that for central AB collisions 
the CERN data require the introduction of final 
state interaction of secondaries. 

Apart from this assumption the model has a 
sound theoretical basis and has great predictive 
power. The latter is due to the fact that strings 
are universal. The momentum distribution and 
fragmentation functions, allowing to compute the 
properties of each individual string, are the same 
in all hadronic and nuclear collisions. Moreover, 
the momentum distribution functions are deter- 
mined from Regge intercepts. For instance, that 
of a valence quark in a pion is given by 

Pl( Xqv )cxx-y 2 (i-x qv y 1/2 • (4) 

The behaviour for x — > is the well known 1 / y/x 
behaviour associated to an ordinary Regge tra- 
jectory of intercept 1/2. Since x qv + Xq v = 1, the 
behaviour for x qv — > 1 is the same as for Xq v — ^ 0, 
and therefore one has the symmetric behaviour 
(jj). In this way the end point behaviour is un- 



Figure 2. Four string component in nucleon- 
nucleus scattering for two inelastic collisions. 

der control. Of course, we can multiply ([|) by a 
smooth function of x. However, in view of the sin- 
gular behaviour at the end points, this will have 
little influence on the numerical results. 

For a valence quark in a proton, one has instead 

^(x) ocaT^l-z) 1 - 5 . (5) 

Here the behaviour at x — > 1 is obtained by 
taking the diquark momentum fraction to zero. 
This behaviour is controlled by an exotic qqqq 
Regge trajectory of intercept close to - 3/2. As 
a consequence of (|^) the valence quark in a pro- 
ton is slow in average whereas the diquark is fast 
(xq v + x qq = 1). Its fragmentation into a lea- 
ding baryon accounts, in this way, for the leading 
particle effect (see Section 2.3). 

In the first applications of DPM there was no 
attempt to compute the fragmentation functions. 
Only their universality was used. For instance in 
studying pA collisions, one can obtain the frag- 
mentation function from a fit to pp data (or the 
pA data for a given A). The model then predicts 
both the A dependence and the energy depen- 
dence. An important progress in the model hap- 
pened when Kaidalov showed || that the frag- 
mentation functions can also be determined to a 
large extent from Regge arguments. In this way 



Figure 3. Feynman x distribution of proton and 
antiproton in pp collisions at 200 GeV/c com- 
pared with the data of ref. [30] at 175 GeV/c. 
Diffraction is not included in the calculation. 



the model becomes very predictive and the spec- 
tra of different hadrons can be computed basi- 
cally without free parameters - except absolute 
normalizations . 

In order to illustrate this point, I show in Figs. 
3 and 4 the leading proton and leading A xf dis- 
tribution in pp collisions (J. The shape of the 
leading proton spectrum is practically fixed the- 
oretically (here we show only the non-diffractive 
part but the diffractive one can also be computed 
in a similar way ||). Moreover, the different 
shapes of the proton and A distributions is en- 
tirely due to the difference between a#(0) = 1/2, 
which appears in the power of 1 — xp for the lea- 
ding proton, and 0^(0) = 0, which controls the 
xf — * 1 behaviour of the A-distribution. Like- 
wise, the sharp decrease of p and A distributions, 
also shown in Figs. 3 and 4, is entirely controlled 
by known Regge intercepts Q , [|| . 

2.2. Nucleus-nucleus collisions 

The generalization of the model to nucleus- 
nucleus collisions is rather straightforward. For 



Table 1 

The number of charged secondaries per unit ra- 
pidity at y* — for average (po" er ) and central 
collision (p^ ntral ) a t RHIC and LHC [6]. 
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Figure 4. Feynman x distribution of A and A pro- 
duced in pp collisions at 200 GeV/c are compared 
to a compilation of data from ref. [5] . 



simplicity I consider here A A collisions. In the 
approximation of only two strings per nucleon- 
nuclcon collision, the rapidity distribution of se- 
condaries is given by 



dN 



AA 



dy 



-(y) = n A 



N 



11 A, 



-lAr, 



(y) + N** 



+ 2{n-n A ) N"(y) 



-11A T{y) 

(0) 



Here N(y) are the rapidity distributions of the 
individual strings, n A is the average number of 
wounded nucleons of A and n is the average num- 
ber of nucleon-nucleon collisions. Both ha and 
n can be computed in the Glauber model. For 
instance for an average collision (i.e. integrated 
over impact parameter), one has 



A 2 a NN /a AB oc A 4/3 



(7) 



respectively) and as many valence quarks. This 
accounts for the first term of (||). The remai- 
ning strings : 2(n — n A ) have to be stretched by 
sea quarks and antiquarks, because the availa- 
ble valence constituents are all included in the 
first term ; this accounts for the second term of 
Of course we should combine the valence 
and sea constituents of the projectile with those 
of the target in all possible ways. However, for 
linear quantities such as multiplicities, each or- 
dering gives practically the same result. 

We can see from (JgJ) and ((?]) that, if all strings 
would have the same plateau height (i.e. the same 
value of iV(0)), the plateau height in an average 
AA collision would increase like A 4 / 3 . However, 
at present energies, the plateau height of the q s - 
q s , strings is smaller than that of qq-q ones, and 
the first term of ([|) dominates. At higher ener- 
gies the contribution of the sea strings becomes 
increasingly important. Therefore, in order to 
make predictions for RHIC and LHC we have to 
introduce the multistring configurations in each 
nucleon-nucleon collision. If their average num- 
ber is 2K (this number can be computed in a 
generalized cikonal model ; one gets K ~ 2 at 
y/s = 200 and K w 3 at yfs = 7 TeV) the total 
number of strings is 2Kn, and eq. (@) is changed 
into 



dN 



AA 



dy 



<y) 



n A 



N iq Ap -q v AT ( y } + N<iA P - 



HAn 



(y) 



Note that the total number of strings is 2n, i.e. 
two strings per inelastic nucleon-nucleon collision. 

The interpretation of (^|) is obvious. With t%a 
struck nucleons, we have at our disposal n A di- 
quarks of projectile and target (qq Ap and qq Ar , 



(2K - 2)N q °-1°] +(n- n A )2KN q " 



(8) 



The predictions || for average and central (b w 
0) SS and Pb-Pb collisions at RHIC and LHC 
energies are given in Table 1. 



Figure 5. Same as Fig. 2 for the diquark breaking 
component. 



2.3. Diquark breaking and the leading 
baryon 

^From Eq. (0) it follows that the diquark is 
fast in average and carries a large longitudinal 
fraction of the incoming proton. Its fragmenta- 
tion will produce a leading baryon as illustrated 
in Figs. 1 and 2. These configurations corres- 
pond to the fragmentation of the diquark as a 
whole (the two valence quarks of the diquark find 
themselves into the same final state baryon) , and 
will be called in what follows diquark preserving 
(DP) components. However, the diquark can also 
break in the way illustrated in Figs. 5 and 6, 
which we shall call in what follows diquark brea- 
king (DB) components. It was stressed in M that 
these two components give contributions to the 
central plateau height which have different ener- 
gy behaviour. Therefore when studying pp col- 
lisions in a broad range of energies, one should 
keep track separately of these two components. 
This is even more necessary in nuclear collisions 
since the two components have very different A- 
dependence (the ratio of DB over the DP compo- 
nents increases with ^4). This is very important 
since the DB component produces baryons which 
are slower than those produced by the DP one. 

In order to simplify the presentation I discuss 



Figure 6. Same as Fig. 5 but here the produced 
baryon is formed by three sea quarks. 



pA collisions. Following ||, I divide the total 
nucleon-nucleon cross-section into a DP piece 
and a DB one : <Ti n = <tdp + udb- I assume 
that once the diquark has been destroyed in a col- 
lision with one nucleon of the nucleus it cannot 
be reconstructed in further collisions with other 
nucleons of the nucleus. The NA cross-section 
involving n inelastic NN collisions is then given 

by 

x [l-a in T A (b)] A ~ n . (9) 

Here Ta(J}) is the standard nuclear profile 
function at impact parameter 6, normalized to 
unity. In Eq. (^J) we have replaced the 
usual factor cr-^TJ corresponding to the cross- 
section for n inelastic collisions by the product 

n 

£ ObB^pT'I = « - P%p)T2(b). Indeed 

i=l 

only the term cr^pT^ will contribute to the di- 
quark preserving cross-section. Summing in n we 



collisions. We have [jlO| 



^%(b) = l-(l-<J DB T AB {b)) 



AB 



(12) 



Figure 7. Rapidity distribution of the p-p diffe- 
rence in peripheral and central SS collisions (solid 
curves) [8]. The dashed curves are obtained with 
(Jdb = 0. The data are from [11]. The dot- 
ted curve is our prediction for the nucleon minus 
antinucleon rapidity distribution in central Pb-Pb 
collisions. 



have 



N A 
'DB 



(b) = £ = Ml - <TD B T A (b)] A -(io) 



n=l 



Eq. (|10|) shows that the diquark preserving cross- 
section belongs to a class of processes which 
has only self- absorption (or self-shadowing). Ob- 
viously 



4Mb) 



NA 



(b)-<7D A B (b) 



(11) 



Since o\db < "in, it is clear from Eqs. (10) and 



(^TJ) that (Jpg increases with A faster than cr B p. 
Actually, when gdb is sufficiently small to neglect 
in Eq. (|l^) second and higher powers of ctdb, 
a DB wm increase linearly with A. This proves 
the result stated above that the relative size of 
the DB component increases with increasing A. 
The result can be easily generalized to an AB 



where T AB (h) = J d 2 sT A (s)T B (b - s). For a DB 
sufficiently small we get from (|l2|), after integra- 
tion in impact parameter, <x^,g = ABa BB . 

In the numerical calculations, I take <j bb = 7 
mb corresponding to a 20 % weight of the DB 
component. Using eq. (|l2|) one finds that this 
weight has increased to 40 % for central SS col- 
lisions and is as large as 64 % in central Pb-Pb 
collisions. 

The results for peripheral and central SS col- 
lisions at CERN energies are shown in Fig. 7 
and compared with the data from the experi- 
ment NA35 PH . We also show the results ob- 
tained with gdb = 0. We see that the dramatic 
dip present in the latter case for central SS colli- 
sion has been largely filled in by the contribution 
of the DB component in agreement with experi- 
ment. The prediction for a central Pb-Pb collision 
is also shown. In this case the dip is converted 
into a broad plateau - in agreement with prelimi- 
nary data from the NA49 collaboration fll^ . 

2.4. DPM and QGP signals 

In the following I will discuss a possible inter- 
pretation of two QGP signals (strangeness en- 
hancement and J ftp suppression) in the frame- 
work of DPM as well as the required modifications 
of the model - such as string fusion and final state 
interaction of secondaries. I will also discuss pos- 
sible consequences of string fusion for cosmic ray 
physics. 

2.5. Strangeness enhancement 

How can one explain in DPM the enhancement 
of the ratios (|J) ? Since the average number of 
strings increases with A (or with centrality), and 
since the extra strings involve sea quarks and an- 
tiquarks at their ends, strangeness will be en- 
hanced provided the fraction of s-quarks in the 
sea is larger than the ratio of strange over non- 
strange particles at present energies. However, 
this mechanism is numerically important only for 
kaons. A/A production in a q s -q s string is negli- 
geable at present energies due to threshold effects. 
A possible way to enhance A/ A production is via 



string fusion. When two strings overlap, both in 
transverse space and in rapidity, their momenta 
as well as their quantum number can merge. If 
for instance two strings us and sd fuse into a us- 
sd one, it will be easy to produce a A/ A pair. 
In this way it has been possible to explain A en- 
hancement [[[3). 

While this is a very plausible mechanism, the 
evidence for string fusion is not yet compelling. 
Indeed, one can assume that diquark-antidiquark 
pairs are present in the proton sea (with the same 
relative amount required in the string breaking 
process in order to produce baryon pairs) [fbif . 
This mechanism is equivalent to string fusion in 
order to explain A enhancement but does not have 
some of the consequences of string fusion that will 
be discussed later. 

A limitation of the string fusion mechanism is 
that it gives the same enhancement of A and A in 
absolute value. Experimentally the enhancement 
of A is much bigger. What can produce it ? First, 
the DB component introduced in Section 2.3 is 
a source of A-enhancement at mid-rapidities. In- 
deed, in the DB diagram of Fig. 6 the probabi- 
lity to produce a A is three times larger than in 
the DP component of Fig. 2 since the strange 
quark can be any of the three sea quarks that 
form the baryon. For the DB diagram of Fig. 
5 this probability is two times larger than that 
of the DP component. However, the calcula- 
tion shows that the obtained A-enhancement is 
substantially smaller than the experimental one 
|l5| . Another source of A enhancement is there- 
fore needed. Such a source is provided by fi- 
nal state interaction of the produced secondaries. 
The most important interactions turn out to be 

HIM 



N -»• K + A(K* + A, • • •) 
N ^ K + A(K*+A, •••) 



(13) 
(14) 



Since the strange particle yield is proportional to 
the product of densities of the interacting parti- 
cles, it is clear that A enhancement will be more 
important than that of A due to the larger density 
of nucleons as compared to antinucleons. More- 
over, the reactions 



Table 2 

A rapidity distribution in central SS collisions at 
200 GeV/c per nucleon computed with and with- 
out final state interaction (FSI) compared to data 
from the NA35 collaboration [17]. 



V* 


No FSI 


With FSI 


NA35 





0.78 


1.9 


2.2 ± 0.3 
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0.77 


1.8 


2.1 ± 0.3 
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0.75 


1.7 


2.1 ± 0.3 


1.5 


0.70 


1.6 


2.2 ± 0.3 
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0.62 


1.3 


1.4 ± 0.3 



Table 3 

A rapidity distribution in central Pb Pb collisions 
at 160 GeV/c per nucleon computed with and 
without final state interaction (FSI). 



No FSI With FSI 






8.4 


23 - 


- 31 


0.5 


8.1 


22 - 


- 30 


1 


6.8 


20 - 


- 25 


1.5 


5.0 


16 - 


- 19 


2 


2.9 


8.8 - 


- 9.1 



which produce a decrease of A and A, are less im- 
portant than (|l^)-((lj) due to the smaller density 
of A (A) as compared to nucleons (antinucleons). 

Calculatio ns alon g these lines have been per- 
formed in 



5 16 



A-> K + N 



A ^ K + N , 



(15) 



The results for central SS 
collisions are given in Table 2 and compared with 
experiment jl7|]. The agreement is reasonably 
good. We see that final state interaction is needed 
in order to reach agreement with experiment. 
Predictions for central Pb Pb collisions are given 
in Table 3. The value with final state interac- 
tion is consistent with preliminary data from the 
NA49 collaboration @ (23 ± 7 at y* ~ 0). 

2.6. J/- suppression 

We turn next to J/ip suppression p9[ . For all 
pA and AB data [^0| involving a light projec- 
tile (up to S) there is an alternative explana- 
tion to Debye color screening proposed in Jl9| . 
This explanation is based on the following me- 
chanism ||2ll| . In the present experimental condi- 
tions the J/ijj is produced outside (behind) the 



Figure 8. The ratio fl MM o'( J/ip)/a(DY) versus 
the interaction length L for pp, pA, SU, and Pb 
Pb interactions. Data are from refs. [20,22]. The 
straight line is obtained with nuclear absorption 
alone. The points labelled theory are obtained 
[25] with nuclear absorption plus final state in- 
teraction of the J/ip bound state. 



nucleus. However, the cc pair is produced inside 
the nucleus and can interact with nucleons in its 
path through the nucleus. This interaction can 
modify its wave function in such a way that it 
has a vanishing projection into the J/ip bound 
state. A single absorptive cross-section, a a b s ~ 
6.2 mb, allows to describe all data except the 
Pb Pb ones. This is illustrated in Fig. 8. The 
data for Pb Pb, also displayed in this figure, 
show clearly that nuclear absorption alone does 
not work. Another "conventional" mechanism is 
the final state interaction of the J / ip bound state 
with other secondaries produced at comparable 
velocity (co-movers) producing open charm : 

TT + J/ip^ D + D + -- (16) 



This idea was introduced in the literature a long 
time ago [ p3| . However, recent theoretical calcu- 
lations have shown that the w + J/ip and N + J/ip 
cross-sections increase very slowly from threshold 
(in contrast to ip' + ir which has a very rapid in- 
crease) pi}] . Therefore in the case of the J/ip 
the co-mover cross-section is expected to be very 
small (not larger than 0.5 mb). Because of that 
the idea of J/ip suppression via interaction with 
co-movers has been progressively abandoned. 

In a recent paper [p5| we have shown that in- 
troducing the final state interaction ( |l6|) with a 
cross-section of 0.4 mb, one can explain the recent 
Pb Pb data. The results are shown in Fig. 8. A 
similar result has been obtained independently in 

II- 

One should note, however, that when introdu- 
cing the interaction with co-movers one has to 
reduce the value of the cross-section for nuclear 
absorption. In this way the J/ip yield in pA colli- 
sions does not fall with increasing A as steeply as 
the data seem to indicate (the "effective" absorp- 
tive cross-section resulting from the combined ef- 
fect of nuclear absorption and interaction with 
co-movers is 5 mb while experimental data re- 
quire 6.2 ± 0.7 mb). Moreover the J/ip yield in 
SU falls somewhat faster than indicated by ex- 
periment. At present, however, the data are not 
precise enough to rule out such an explanation 
and further study is needed to distinguish it from 
a QGP scenario |2lJ and/or other collective ef- 
fects such as string percolation |2Sj. 

There are, of course, many other interesting re- 
sults obtained at CERN which I have not men- 
tioned. In particular Bose-Einstein correlation 
measurements (interferometry) indicate that the 
size of the interaction region is larger than the 
geometrical size of the colliding system. Again, 
this requires the existence of final state interac- 
tion of secondaries. 

2.7. Consequence of string fusion for cos- 
mic ray physics 

As discussed in section 2.5, there are hints in 
the data in favor of string fusion, although there 
is no compelling evidence yet. Nevertheless, it is 
interesting to examine its possible consequences 
in cosmic ray physics. In a recent paper the San- 



that these two effects can change the profile of 
extensive air showers and in this way alter con- 
clusions on the primary composition and energy. 
For example, provided their composition could be 
determined by an independent method, cosmic 
ray events with energies around 3.10 20 eV would 
contain, in the presence of string fusion, parti- 
cles with \xf\ > 2 or 3 (or more in the presence 
of string percolation ) , and would thus corres- 
pond to an energy for the primary upto 2~A times 
smaller than the one without string fusion. Such 
a reduction in energy could make these events 
compatible with the cut-off due to the scattering 
of cosmic rays with the microwave background. 



Figure 9. Total average multiplicity versus pri- 
mary energy for a fixed composition assuming 
string fusion (solid line) and with no string fu- 
sion and a composition that changes with energy 
(dashed line). See text for details. 



tiago group |29J discusses two such consequences : 
the rise of the average shower depth of maximum 
for cosmic rays in the energy range 10 16 to 10 19 
eV and the possibility that string fusion acts as a 
hadronic accelerator. 

When the energy increases, string fusion be- 
comes more important and produces a decrease 
of the average multiplicity. Taking a fixed chemi- 
cal composition (i.e. independent of energy), con- 
sisting of 90 % of iron and 10 % of protons, one 
gets an average multiplicity, which, in the range 
10 17 to 10 19 eV, is identical to the one obtained 
without string fusion and a chemical composition 
of primaries changing with energy (from 75 % Fe 
and 25 % p at 10 16 eV to 50 % Fe and 50 % p at 
10 19 eV). This is shown in Fig. 9. 

Another interesting consequence of string fu- 
sion is that it produces secondaries outside the 
kincmatical range of a nucleon-nucleon interac- 
tion. The authors can explain in this way the 
so-called cumulative effect. Moreover, they claim 



3. CONCLUSIONS 

The QGP search has gained full credibility with 
the lattice QCD results. Experimental data are 
of very high quality but give only some hints of 
QGP formation. Future data and especially fu- 
ture machines will be decisive in this search. 

In the mean time, we have gained new insight 
in the dynamics of multiparticle production al- 
lowing to improve current models (strong baryon 
stopping hinting at diquark breaking, string fu- 
sion with interesting consequences for cosmic ray 
physics and final state interaction of secondaries). 

Finally, I want to stress once more that the 
CERN heavy ion program has produced very 
many interesting results not mentioned here. 
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